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We compare the tau neutrino flux arising from the galaxy and the earth atmosphere for 10® < 

E/GeV < The intrinsic and oscillated tau neutrino fluxes from both sources are calculated. 

The intrinsic galactic flux {E > 10® GeV) is calculated by considering the interactions of high- 
energy cosmic-rays with the matter present in our galaxy, whereas the oscillated galactic Ur flux 
is coming from the oscillation of the galactic flux. For the intrinsic atmospheric Vt flux, we 
extend the validity of a previous calculation from E < 10® GeV up to A < 10®*^ GeV. The oscillated 
atmospheric Ur flux is, on the other hand, rather suppressed. We find that, for 10® < E/GeV < 5 • 

10^, the oscillated Ut flux along the galactic plane dominates over the maximal intrinsic atmospheric 
i^T flux, i.e., the flux along the horizontal direction. We also briefly mention the presently envisaged 
prospects for observing these high-energy tau neutrinos. 

PACS numbers: 95.85.Ry, 13.85.Tp, 98.38.-j, 14.60.Pq 


I. INTRODUCTION 

Searching for high-energy tau neutrinos {E > 10® GeV) will yield quite useful information about the highest energy 
phenomenon occurring in the universe Q. The same search may also provide evidence for physics beyond the standard 
model 1^]. The latter is suggested by the recent measurements of the atmospheric muon neutrino deficit, though yet 
there is no observation of oscillated atmospheric tau neutrinos at a significant confidence level . Interestingly, it is 
also only recently that we have the first evidence of existence of tau neutrinos . 

The high-energy tau neutrinos can be produced in pp and pj interactions taking place in cosmos. These interactions 
produce unstable hadrons that decay into tau neutrinos. In this paper, we mainly concentrate on pp interactions and 
will only briefly comment onpy interactions as source interactions for producing high-energy tau neutrinos. There can 
be several astrophysical sites where the pp interactions may occur. Examples of these include the relatively nearby and 
better known astrophysical sites such as our galaxy and the earth atmosphere, where the basic pp interactions occur as 
pA interactions. The pp interactions in these sites form a rather certain background to the extra-galactic high-energy 
tau neutrino searches. It is possible that such interactions are the only sources of high-energy tau neutrinos should 
the search for high-energy tau neutrinos originating from several proposed distant sites such as AGNs, GRBs, as 
well as groups and clusters of galaxies, turns out to be negative. Therefore, it is rather essential to investigate the 
high-energy tau neutrino flux expected from our galaxy and the earth atmosphere. 

Both the galactic and atmospheric tau neutrinos can be categorized into intrinsic and oscillated ones. Here, intrinsic 
Vr flux refers to the Vr produced directly by an interaction while oscillated Vt refers to the resulted from the 
oscillation. Presently, there exists no estimate for the intrinsic high-energy flux originating from our galaxy in pp 
interactions, although the estimates for zzg and fluxes due to the same interactions are available Q . In this work, 
we calculate the intrinsic z/,- flux from our galaxy by using the perturbative and nonperturbative QGD approaches 
to model the pp interactions, and taking into account all major tau neutrino production channels up to A < 10®® 
GeV. We note that the production of tau neutrinos in the terrestrial context was discussed in Ref. |^, which uses a 
nonperturbative QGD approach for pp interactions. To calculate the oscillated galactic and atmospheric Vr fluxes, we 
apply the two-flavor neutrino oscillation analysis 

It is essential to compare the above galactic tau neutrino flux with the flux of atmospheric tau neutrinos. The 
intrinsic atmospheric tau neutrino flux has been calculated for E < 10® GeV Q. In this work, we extend the 
calculation up to A < 10®® GeV. Such an extension requires the input of cosmic-ray flux spectrum for an energy range 
beyond that considered in Ref. 0. Furthermore, for a greater neutrino energy, the solutions of cascade equations 
relevant to the neutrino production behave differently. For the oscillated z/,- flux, it is interesting to note that the 
oscillation length for zz^ ^ Vr for the energy range 10® < E/GeV < 10®® is much greater than the thickness of the 
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earth atmosphere. Hence the oscillated atmospheric flux in this case is highly suppressed. 

One may argue that the interaction of the high-energy cosmic-rays with the ubiquitous cosmic microwave background 
(CMB) photons present in cosmos {pj rather than pp interactions) could also be an important source for high-energy 
astrophysical tau neutrinos. The center-of-mass energy (-^ 5 ) needed to produce a r lepton and a Vr is at least ^ 1.8 
GeV. In a collision between a proton with an energy Ep and a CMB photon with an energy , the center-of-mass 

energy squared of the system satisfies < s < AEpE^^y^^ + nip. Since the peak of the CMB photon flux spectrum 
with a temperature ~ 2.7 K is at about 2.3 • 10“"* eV, it requires a very energetic proton with Ep > 2.5 • 10^^ GeV in 
order to produce a tVt pair. Thus, the contribution of the intrinsic tau neutrino flux from the interaction between the 
cosmic proton and the CMB photon is negligible, unless we are considering extremely high-energy protons, beyond 
the presently observed highest energy cosmic-rays |^. To compute the oscillated Vr flux in this case, we consider 
the non-tau neutrino flux generated by the p"f interaction via the A resonance, commonly referred to as the Greisen- 
Zatsepin-Kuzmin (GZK) cutoff interaction, which assumes that the proton travels a cosmological distance |P. A recent 
calculation of the intrinsic non-tau GZK neutrino flux indicates that this flux peaks typically at A ^ 10® GeV 
beyond the reach of presently operating high-energy neutrino telescopes such as AMANDA and Baikal experiments 
[Q. This flux decreases for E < 10® GeV. In fact, it falls below the intrinsic non-tau galactic-plane neutrino flux for 
E < 5 - 10^ GeV. The neutrino flavor oscillations of the intrinsic non-tau GZK neutrinos into tau neutrinos result in a 
Vr flux comparable to the original non-tau neutrino flux ||ll| . Therefore, in the absence (or smallness) of tau neutrino 
flux from other possible extra-galactic astrophysical sites, the only source of high-energy tau neutrinos besides the 
atmospheric background is from our (plane of) galaxy, typically for 10® < E/GeV < 5 • 10^. This is an energy range 
to be explored by the above high-energy neutrino telescopes in the near future. 

The organization of the paper is as follows. In Section II, we discuss the calculation of intrinsic high-energy tau 
neutrino flux from our galaxy, including the description of the flux formula, the galaxy model, and the various tau 
neutrino production channels taken into account in the calculation. Although this flux will be shown small, we shall 
go through some details of the calculation since they are also useful for the calculation in the next section. In Section 
III, we present our result on the intrinsic atmospheric Vt flux and compare it with the galactic one. In Section IV, we 
discuss the effects of neutrino flavor mixing, which are used to construct the oscillated Vr fluxes from the galaxy and 
the earth atmosphere respectively. The total galactic Vt flux (the sum of intrinsic and oscillated fluxes) is compared 
to its atmospheric counterpart, and the dominant energy range for the former flux is identified. We also mention the 
currently envisaged prospects for identifying the high-energy tau neutrinos. We summarize in Section V. 


II. THE INTRINSIC GALACTIC TAU NEUTRINO FLUX 


A. The tau-neutrino flux formula and the galaxy model 


We use the following formula for computing the Vr flux: 


dV., 

dE 



dEp (j)p{Ep) f{Ep) 


dnpp^iy^jf-Y 

dA 


In the above equation, E is the tau neutrino energy and the cosmic-ray flux spectrum, (j)p{Ep), is given by 


r I.7(Ap/GeV)-®-^ iovEpKEo, 
\ 174 (Ap/GeV)-® for Ep > Eq, 


( 1 ) 

( 2 ) 


where Eq = 5 - 10® GeV and (j)p{Ep) is in units of cm“® s“® sr“® GeV“^. We assume directional isotropy in (j)p{Ep) for 
the above energy range. The function f{Ep) is equal to R/Xpp{Ep), where \pp{Ep) = (CTp“*np)“^ is the pp interaction 
length and i? is a representative distance in the galaxy along the galactic plane. The target particles are taken to 
be protons with a constant number density of 1 cm“® and R is taken to be ^ 10 kpc, where 1 pc ~ 3 • 10^® cm. 
The is the total inelastic pp cross section. Since the high-energy protons traverse a distance R much shorter 
than the proton interaction length, the proton flux spectrum (j)p{Ep) is assumed to be constant over the distance 
R. Furthermore, we calculate the intrinsic tau neutrino flux along the galactic plane only to obtain the maximal 
expected tau neutrino flux. The matter density decreases exponentially in the direction orthogonal to the galactic 
plane, therefore the amount of intrinsic tau neutrino flux decreases by approximately two orders of magnitude for 
the energy range of our interest. For further details, see Ingelman and Thunman in Ref. Q. The function f{Ep) in 
Eq. (|^ basically gives the number of inelastic pp collisions in the distance i?, while the distribution dn/dE is the 
differential cross section normalized by i.e., 

dUpp_,i/^_l_y _ 1 d(Jpp_,i/^_|_y 

dE “ dE 


( 3 ) 
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The above distribution gives the fraction of inelastic pp interactions that goes into z/t-’s. We can now simplify Eq. (|^) 
into 

= Rn, 1 “ dE, d>,(E,) , ( 4 , 

It is clear that the task of computing dN^^/dE relies on the evaluation of the differential cross section da/dE in 
pp interactions. In this work, we include all major production channels of tau neutrinos, namely, via the meson, 
6 -hadron, tt, W* , and Z*. In general, the symbol i^r shall be used to account for the contribution of and Dt unless 
otherwise mentioned. We note that the heavy intermediate states, such as the meson, 6-hadrons and other heavier 
states, decay (into i^,-) before they interact with other particles in the galactic plane. This is due to the rather small 
matter density of the medium and the large distance between the proton source and the earth. Before we proceed, 
let us remark that the intrinsic tau neutrino production by the galactic py interactions is suppressed relative to that 
in the galactic pp interaction, because <C Up for the energy range of interest to us p^. 


B. Tau neutrino production 

1. Via Ds mesons 

The lightest meson that can decay into a t-Vt pair is the Ds meson. It was pointed out that the production and 
the decay of Dg meson is the major production channel for tau neutrinos in the AGN 1^. We expect the same will 
be true for galactic tau neutrinos. The Ds meson decays into a charged r lepton and a Vt- The charged t lepton 
subsequently decays into the second Vr plus other particles, which can be one prong or three prong. The kinematics 
of the T lepton decay is treated by the Monte Carlo technique. For simplicity, we assume that the r lepton decays 
into a Vr and a particle Y with the mass my satisfying 0.1 GeV < my < rrir — 0.1 GeV. Consequently, the second Vr 
is much more energetic than the first one because the Dg mass is only slightly larger than rur ■ We take the branching 
ratio B{Ds —> as ~ 0.074 |l^ . 

We now turn to the production of Dg mesons. Here, we employ two approaches to calculate the production of Dg 
mesons: (i) the perturbative QCD (PQCD) and (ii) the quark-gluon string model (QGSM). In the PQCD approach, 
we use the leading-order result for pp cc: 

a{pp ^ cc) = Y^ j J da:ida ;2 fi/p{xi)fj/p{x 2 )d{ij cc) , (5) 

where fi/p{x) are the parton distribution functions (we use the CTEQv5 [p^ ), while the parton subprocesses are 
cc. We use a K factor, AT = 2, to account for the NLO corrections |Q. The matrix elements for these 
subprocesses can be found in Appendix A. The c or c then undergoes fragmentation into the Dg meson, which we 
model by the Peterson fragmentation function ||^ with e « 0.029 (see Appendix A). The probability of a 

charm quark fragmenting into a Dg is 0.19 |Q (we have added the fc^D^ and fc—>D* )- Typically, a fraction z (z < 1) 
of the charm-quark energy is transferred to the Dg meson so that the energy spectrum of Dg is softer than that of 
the charm quark. 

The QGSM approach is nonperturbative and is based on the string fragmentation. It contains a number of parame¬ 
ters determined by experiments . The production cross section of the Dg meson is given by the sum of n-pomeron 
terms 


dcr^'’(s,a;) 

dx 






( 6 ) 


where x = 2p\\/^/s and xy = -\-p\)/s. The functions and (j)^‘‘{s,x) are given in Appendix B. 

A comparison of these two approaches for Dg meson is shown in Fig. |^. The dashed line in the figure is the spectrum 
of the injected proton flux given by Eq. (^). The Vr spectra calculated by these two approaches agree well with each 
other for E < 10® GeV. Beyond this energy, the QGSM approach gives a relatively harder spectrum. In fact, this 
behavior was already seen in the da/dx distribution. Nevertheless, in the region where the two approaches differ, 
the tau neutrino flux is already small. To our knowledge, the heaviest meson production that the QGSM has been 
applied to is the Dg meson production. The current highest energy collider experiment for Dg meson production is at 
the FERMILAB TEVATRON with d/s = 1.8 • 10® GeV, which corresponds to an Ep ~ 1.7 • 10® GeV, as s ~ 2mpEp 
in our setting. Note that up to this d/s, the agreement between the two approaches is quite good, according to Fig. 
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FIG. 1: A comparison between the PQCD and QGSM approach to the energy spectrum of the intrinsic galactic Vt flux coming 
from the Da meson. The thick dashed curve is the injected proton flux spectrum given by Eq. 


1^. We have used a factorization scale = s/4 and the one-loop running strong coupling constant as with the value 
ctsiQ^ = = 0.118, and the rUc = 1.35 GeV. 

An important quantity in the neutrino flux calculation is the average fraction of the injected proton energy being 
transferred to the tau neutrino, i.e., the ratio r = E/Ep. The average value of r is given either by the mean 



( 7 ) 


or by the value of r at which the distribution da/dE attains the peak. We have found that both averages of r are 
very close to each other. The (r) ranges from 5 • 10“^ to 5 • 10“^ for Ep from 10^ to 10^^ GeV for the production 
channel pp ^ cc ^ Dg + Y ^ + Y, using the PQCD approach. The higher the injected proton energy, the smaller 

is the fraction of the incident Ep that goes into hadrons. 


2. Via bb, tt, W* and Z* 

The production of bb and tt in pp interactions can be calculated quite reliably by the PQCD approach, similar to 
the calculation of cc. The relevant matrix elements, including the ones for W* and Z*, are listed in Appendix A. The 
results ^are shown in Fig. |2[ A few observations can be drawn from the figure pl| . (i) The production via Dg mesons 
dominates for E < 10® GeV, followed by 6-hadrons, W*, Z*, and tt respectively, (ii) For E > 10® GeV all these 


^ After taking into account the sources of uncertainties such as values of rric, nif,, nit, R, np, etc., we estimate that our calculation of the 
intrinsic galactic Ut flux is reliable within an order of magnitude. 
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FIG. 2: Intrinsic galactic tan neutrino flux calculated via various intermediate states and channels: Da, 6-hadron, W*, Z*, 
and tt. The injected proton flux spectrum is also shown. 


production channels become comparable, (hi) The intrinsic tau neutrino flux is about 10 — 12 orders of magnitude 
smaller than the injected proton flux. 


III. THE INTRINSIC ATMOSPHERIC TAU NEUTRINO FLUX 


The earth atmosphere is an interesting extra-terrestrial site where the basic pp interaction occurs in the form of a 
pA collision with A the nuclei present in the earth atmosphere. Incidently, it is the only known nearby extra-terrestrial 
site from where the intrinsic neutrinos are observed as a result of high-energy cosmic-ray interactions. 

We have calculated the downward and horizontal intrinsic atmospheric flux for the energy range 10^ < E/GeV < 
10^^. We used the nonperturbative QCD approach mentioned in the last section to model the production of Dg mesons 
in pA interactions. We have used the (j)p{Ep) given by Eq. (||) and the Z-moment description for the calculation of 
intrinsic tau neutrino flux, which is appropriate for a varying target density medium . Since the tau neutrino flux 
is determined by the flux of Dg meson, we briefly discuss the cascade equation for the Dg flux. In general, we have 


<^4>Da 

dX 


4 >Ds 

^Ds 


'pDs 

pdos 


I 7 ^De I 7 

+ ZDsDs T-r ZpD,—, 

^Ds 


( 8 ) 


where X is the slant depth, i.e., the amount of atmosphere (in g/cm^) traversed by the Dg meson {X = 0 at the top 
of the atmosphere), (j)p{E,X) and (j)Da{E,X) are the fluxes of protons and Dg mesons respectively. The Ad„ and doa 
are the interaction thickness (in g/cm^) and the decay length of the Dg meson respectively. Finally, the Z-moments 
ZpDa s-nd describe the effectiveness of generating Dg meson from the higher-energy protons and Dg mesons 

respectively. We have 


ZpDa{E) 

ZOaDaiE) 



.j:,/ (t>p{E',0) Ap(E) dnpA^Da+Y{E,E') 

ME,0) Ap(E') dE 

,j:,i 4’DsiE',0) XpajE) dnp^A^Ds+YiE, E') 
(l>Da{E,0) XoaiE') dE 


E 


(9) 
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where dnp^_>_D^+y/di? and Aud^a^d^+y/<^E are defined according to Eq. (^). We note that Eq. (||) should be solved 
together with the cascade equation governing the propagation of high-energy cosmic-ray protons. In fact, the proton 
flux equation can be easily solved such that 

« exp(-^)(/)p(£;,0), (10) 

l\p 

where Ap = Ap/(1 — Zpp) is the proton attenuation length with Ap the proton interaction thickness (in g/cw?) and 
the Zpp given by 



XpjE) dnpA^p+YiE,E') 
(^p(E;,0) Ap(A0 d£; 


( 11 ) 


An analytic solution of Eq. (^ can be obtained for either the low or the high energy limit. Such limits are characterized 
by whether the decays before it interacts with the medium or vice versa. The critical energy separating the two 
limits is approximately ~ 8.5 • 10^ GeV. In the low energy limit, we disregard the first and third terms in the R.H.S. 
of Eq. (H). On the other hand, one can drop the second term in the high energy limit. With (j)u^ determined, the i^r 
flux can be calculated by considering the decay Dg ^ VtT and the subsequent decay t ^ Vr + Y as treated in Section 
II. We first obtain two i^r fluxes, valid for low and high energy limits respectively, in terms of Z-moments and then 
interpolate the two fluxes. A complete numerical solution without the interpolation is given in a separate work 



Log,„(E/GeV) 


FIG. 3: Intrinsic horizontal Ur flux via production and decay of the Ds meson in the earth atmosphere. For 10® < E /GeV < 10®, 
the results by PR are also shown. In the inset, we compare our calculated ZpDs with the one given by rescaling TIG’s Zpjjo. 
The total intrinsic galactic-plane tau neutrino flux is also shown. 

In Fig. 1^, we show our result for the intrinsic atmospheric flux along the horizontal direction. For comparison, 
the results by Pasquali and Reno (PR) H, valid for 10® < E/GeV < 10®, are also shown. The Vp flux along all 
the other direction is small. For example”, the downward flux is about 8 times smaller than the horizontal one 


^ The upward going high-energy Up with a energy E > 10'* GeV is degraded in energy to about E < 10® GeV after crossing the earth. 
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for E > 10® GeV. We remark that the major uncertainty for determining the above Vr flux is the Z-moment . 
In Ref. 0, the authors calculate using two different approaches, which then give rise to different results for 

the Vr flux. The first approach is based upon next-to-leading order (NLO) perturbative QCD |Q, while the second 
approach rescales the .Zp^o given by the PYTHIA j2^ calculation of Thunman, Ingelman, and Gondolo (TIG) [^ . 
In the inset of Fig. ^ , we also show our calculated Zpu, in comparison with the one given by rescaling TGI’s result 
for Zp^o. We do not show the ZpD^ obtained by NLO perturbative QCD since it is not explicitly given in Ref. . 


IV. EFFECTS OF OSCILLATIONS AND PROSPECTS FOR OBSERVATIONS 

In the context of two neutrino flavors, Vp and Vr, the total Vr flux, dV);°*/d(logiQ E), is given by 

dA*f/d(logio E) = P- dW^/d(logio E) + {l-P)- diV.Jd(logio E). (12) 

Here P = P{vp Vr) = sin^ 29 ■ sin^(Z/Zosc)- The neutrino flavor oscillation length for Vp Vr is Gsc {E/6m?). 
For 10® < E/GeV < 10^® and with Sm^ ~ 10“® eV^, we obtain 10“® < Gsc/pc < 1. We assume maximal flavor 
mixing between Vp and Vr- 

For intrinsic neutrinos produced along the galactic plane, we take dNi, /d(log]^Q E) given by Ingelman and Thunman 
in Ref. [D by extrapolating it up to A < 10^® GeV, whereas for dV^^/d(logio A), we use our results obtained in 
Section II. For galactic-plane neutrinos, we note that Gsc 'C where I ~ 5 kpc is the typical average distance the 
intrinsic high-energy muon neutrinos traverse after being produced in our galaxy. Eq. (^) then implies that, on the 
average, half of the muon neutrino flux will be oscillated into tau neutrino flux, reducing its intrinsic level to one half. 

For downward going neutrinos produced in the earth atmosphere, we take I ~ 20 km as an example. We use the 
(prompt) dAi,^/d(logiQ E) given in [^, whereas for dA;,^/d(logiQ E), we use our results obtained in Section III. For 
horizontal and upward going atmospheric neutrinos, I ^ 10® — 10“^ km. Here, Gsc ^ I, and so the intrinsic atmospheric 
tau neutrino flux dominates over the oscillated one for E > 10® GeV, essentially irrespective of the incident direction. 
We present these results in Fig. along with the GZK oscillated tau neutrino flux briefly mentioned in Section I. 
For the GZK neutrinos, I > Mpc and dVj,^/d(log]^Q A) is taken from Ref. |^. From the figure, we note that the 
galactic-plane oscillated Vr flux dominates over the intrinsic atmospheric Vr flux for E < 5 ■ 10^ GeV, whereas the 
GZK oscillated tau neutrino flux dominates for E > 5 • 10^ GeV. 

A prospective search for high-energy tau neutrinos can be done by appropriately utilizing the characteristic r lepton 
range in deep inelastic (charged current) tau-neutrino-nucleon interactions, in addition to the attempt of observing the 
associated showers. For E close to 6 TO® GeV, the (anti electron) neutrino-electron resonant scattering is also available 
to the search for (secondary) high-energy tau neutrinos |2^ . The main advantages of using the latter channel are that 
the neutrino flavor in the initial state is least affected by neutrino flavor oscillations and that this cross section is free 
from theoretical uncertainties ||^ . The appropriate utilization of the characteristic r lepton range in both interaction 
channels not only helps to identify the incident neutrino flavor but also helps to bracket the incident neutrino energy 
as well, at least in principle. 

For downward going or near horizontal high-energy tau neutrinos, the deep inelastic neutrino-nucleon scattering, 
occurring near or inside the detector, produces two (hadronic) showers ]M| |. The first shower is due to a charged- 
current neutrino-nucleon deep inelastic scattering, whereas the second shower is due to the (hadronic) decay of the 
associated r lepton produced in the first shower. It might be possible for the proposed large neutrino telescopes 
such as IceGube to constrain the two showers simultaneously for 10® < E/GeV < 10^, depending on the achievable 
shower separation capabilities (see, also, Q). Here, the two showers develop mainly in ice. Using the same shower 
separation criteria as given in ||32|], we note that the ^ (100 m)® proposed neutrino detector, commonly called the 
megaton detector Q, may constrain the two showers separated by > 10 m, typically for 5 • 10® < E/GeV < 10®. 
The two nearly horizontal showers may also possibly be contained in a large surface area detector array such as Pierre 
Auger, typically for 5 • 10® < E/GeV < 10® jsi). In contrast to previous situations, here the two showers develop 
mainly in air. Several different suggestions have recently been made to measure only one shower, which is due to 
the r lepton decay, typically for 10® < E/GeV < 10^®, while the first shower is considered to be mainly absorbed in 
the earth Q. The upward going high-energy tau neutrinos with E > 10^ GeV, on the other hand, may avoid earth 
shadowing to a certain extent because of the characteristic r lepton range, unlike the upward going electron and muon 
neutrinos, and may appear as a rather small pile up of vi {I = e, andr) with E ^ 10® GeV |^, However, the 
empirical determination of incident tau neutrino energy seems rather challenging here. 

The above studies indicate that for a rather large range of tau neutrino energy, a prospective search may be carried 
out. The event rate in each experimental configuration is directly proportional to the incident Vr flux and the effective 
area of the detector concerned. Presently, no direct empirical upper bounds (or observations) for high-energy tau 
neutrinos exist. 



FIG. 4: Galactic-plane, horizontal atmospheric and GZK tau neutrino fluxes under the assumption of neutrino flavor oscilla¬ 
tions. 


To have an idea of the event rate, let us consider the downward going high-energy tau neutrinos originating from 
the galactic-plane due to neutrino flavor oscillations. The flux of such neutrinos has been given in Fig. 4. Following 
Ref. Q , we note that these galactic tau neutrinos give a representative event rate of < 1 per year per steradian for 
two separable and contained showers with E ~ 10® GeV in a km® volume neutrino telescope such as the proposed 
IceCube. 


V. DISCUSSION AND CONCLUSIONS 

We have calculated the Vt flux due to pp interactions in our galaxy. This flux consists of intrinsic tau neutrino flux 
and that arising from the oscillations of muon neutrinos. We note that the latter flux is dominant over the former by 
four to five orders of magnitude for the considered neutrino energy range. From Fig. 4, one can see that the main 
background for the search of high-energy extra-galactic tau neutrino is due to the muon neutrinos produced in the 
galactic-plane, which then oscillate into tau neutrinos. Such a flux dominates for 10® < E/GeV < 5 ■ 10^. Therefore 
it is clear that searching for extra-galactic tau neutrinos orthogonal to the galactic-plane is more prospective. 

In the calculation of galactic tau neutrino flux, we have used a simplified model of matter distribution along our 
galactic plane to obtain the maximal intrinsic tau neutrino flux. We have explicitly calculated the contribution of 
heavier states such as bb, tt as well as W* and Z* in addition to the more conventional Dg channel to the intrinsic tau 
neutrino flux. We have estimated the average fraction of the incident cosmic-ray energy that goes into tau neutrinos 
and found it to be less than 1%. The contributions from bb, tt, W* and Z* channels are comparable to Dg for E > 10® 
GeV. For Dg channel, we have used both perturbative and nonperturbative QGD approaches. 

We have extended a previous calculation of intrinsic atmospheric lyr flux from E < 10® GeV up to Fi < 10®® GeV. 
Here, we used the nonperturbative QGD approach to calculate the production of Dg mesons in pA interactions. In 
comparison with the intrinsic galactic-plane Vt flux, it is large. However, since the distance between the detector and 
the neutrino source in the galactic plane is sufficiently large, the neutrino flavor oscillations of non-tau neutrinos into 
tau neutrinos makes the eventual tau neutrino flux along the galactic plane greater than the atmospheric tau neutrino 
flux for 10® < E/GeV < 5 • 10®". However, the intrinsic atmospheric flux dominates over the oscillated galactic 
Vr flux in the direction orthogonal to the galactic plane. We have also briefly mentioned the presently envisaged 
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prospects for observations. In summary, we have completed the compilation of all definite sources of tau neutrino 
flux, i.e., those from our galaxy and from the earth atmosphere. Such a compilation is needed before one conducts 
the search for tau neutrinos from extra-galactic sources. 
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APPENDIX A: FORMULAS FOR PQCD 


In this appendix, we list the matrix element squared for the subprocesses of pp —*■ QQ, where Q = c,b, t, used in 
the PQCD calculation. 




d cos 9* 
da 

dcos0* 


727rs3 


(gg^QQ) = 


(m| - if + {m^Q - uf + 2sm\ 
4 


-r —(—rrin — 2>mnt — rrinU + ut) 

Q Q Q ’ 


— - - irr?QU - m%i + ut) 

rriQ) 


{u - m%){t - rrff 


9 1,4 


(2mQ + i+u) + 18^(toq - mqf + u) + iu) 


t — rriQ s 


— {mQ — 2mQi + ui) + 


s 

9 1 


u — rriQ s 


-{rriQ - 2mQU + ut) }, 


(Al) 


where = dTra^, jS = — (4toq/ s) and s, i, u are the usual Mandelstem variables with i—rrf^ = — (s/2)(l—/3cos0*). 

The subprocesses for high-energy tau neutrino production via W* and Z* are qq' —> W* —> and qq ^ Z* ^ 

VtVt- The spin- and color-averaged amplitude squared for these subprocesses are given by 


\M{ud^W* 

Y\M{qq^Z* 



9^ 

1 

i{i — 

12 (s 

-m^)2-pr^TO^ 


9"^ 

1 

u{u - 

12 (s 

- -bP^TO^ 

l^T ^r) 1 

9^ 

1 


3cos^ 

6»w (s - to |)2 -b r| 

■ml 


'nl), 

mf), 

i9l9lfu^ + i9l9'kfP 


(A2) 


f Q 

where g)^ = T^f — Qf sin 9w and 9w is the weak mixing angle, T^f is the third component of the weak isospin and Qf 
is the electric charge in units of proton charge of the fermion /. 

The Peterson fragmentation function is given by 


Dq^Hq^z) = N- 


fl-zf 


[(l-z)2 + ez]2’ 

where N is the normalization constant and e is given in the text. 


(A3) 


APPENDIX B: FORMULAS FOR QGSM 
In Eq. (1^, the functions al^{s) and (j)f‘{s,x) are given as follows: 


4>f‘{s,x) = [Ff^fx+,n)Flf‘{x-,n) + Flffx+,n)Ffq‘{x-,n) + 2{n- l)Ff^^^{x+,n)Ffi‘fx-,n)\ , 


Fffx,n) = ^J dxi f^{xi,n) Gf-’ ^ / dxi/^(a:i, n) , 

Pqq‘i.x,n) = dxi fp‘‘{xi,n) + ^J da;i /““(xi,n) Gf‘ , 


(Bl) 


(B2) 

(B3) 


where 
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= i + 2S. + 25^ {l 

+ / dxi /p"'‘“(a;i,n) 

J X 

+(5s f dxi /p'’''“(a:i,n) 

J X 

+6c f dxi 

J X 


G^‘ ( — 1 + G,rM — 


D, I ^ 

Xi 




D, / X 


Xi 


G^M-l+GfM- 

Xi 


G 




(B4) 


In the above, /p(x,n)’s are the distribution functions describing the n—Pomeron distribution functions of quarks 
or diquarks {i = u, d, uu...) with a fraction of energy x from the proton, and G^( 2 ;)’s are the fragmentation functions 
of the quark or diquark chain into a hadron h which carries a fraction z of its energy. 

The list of the fp{x, n) are given by 

_ r(l+n— 2 q;^) -aR ('1 \aH- 2 a]v + (n-l) 

" r(l-a«)r(a«-2a^ + n)^^ 

_ f«sea 

J P 1 

r(2 + n — 2aN) 

J ' 


fpix,n) = 


r(l - Ur) T{aR - 2ajv + n + 1 ) 

= 

J p 5 


y _ j.'^OLR—2oiN+n 


frM = 
Cm = 


r(2 + n — 2aN) 


^ ’ r(—Ofl + n) r(Q;fl — 2 Q;Ar + 1 ) 


r(l + n — 2aM) 


r(—Ofl + n) r(Q;fl — 2aN + 2) 


X X 


X X 


ajt —2ajv + l ^ 2 ^ _ «« + ("—1) 


(1-x) 


-aR + {n-l) 


fSaealp _ r(l + n + 2aj^ 2 aN 2 a^) /i _ \ 2 aR- 2 aN + {ri-l)-a4. 

^ r(l-a^)r(2a«-2a^ + n-a^) 

fc,ea(p pi) — r(l + n + 2Qff; — 2aN — 2aip) -a^ /1 _ ) 2 aR- 2 aN + {n-l)-a^ 

^ ~ T{l-a^)T{2aR-2aN + n-a^) ^ ’ 


where P is the usual Gamma function. 
The list of the G^ ‘ (z) are given by 


G 


.D = 

u,u,d,d 

^Di 


(z) = 

_ fi \A —a^+ai? —2Q:iv —0:^+2 


GM^) = (1-^) 

= G?»”(z), 


Gf» (z) = (1-zM^ X (1 + ai z"), 

= G^^z), 


G^i(z) = zi-“’^(l-z)^-“^ 
In the above, the input parameters are as follows: 


an = 0.5 , 
ttAT = —0.5 , 

02(ji — 0 , 

otif, = —2.18 , 

A = 0.5 , 

(5^ = 0.25 , 
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= 0.0007 , 
ai = 5 , 

i5c = 0 (0.01), if charm sea contribution is turned off (on). 
The function crPP(s) is given by the following formulas: 

n—1 ^ 

l-exp(-0)^ p- 



’(?) = 


Up 

n z 


(B5) 


where 


z 

(Jp 


In 


1 (GeV)2 


i?2 + Op' ^ 

87r7p exp(^A) . 


exp(^A), 


The best fit parameters are as follows : 
(i) for ^/s < 10^ GeV 


7p 

= 3.64 

(GeV) 


= 3.56 

(GeV) 

Op 

= 0.25 

(GeV)-2 

C 

= 1.5, 


A 

= 0.07. 



(ii) for ^/s > 10^ GeV 


7p 

= 1.77 

(GeV)-2 


= 3.18 

(GeV) 

a'p 

= 0.25 

(GeV) 

C 

= 1.5, 


A 

= 0.139 




